Depinning exponents of the driven long-range elastic string 
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C ' We perform a high-precision calculation of the critical exponents for the long-range elastic string 

driven through quenched disorder at the depinning transition, at zero temperature. Large-scale 
simulations avoid finite-size effects and improve accuracy. We explicitly demonstrate the equivalence 
of fixed-velocity and fixed-driving-force simulations. The roughness, growth, and velocity exponents 
are calculated independently, and the dynamic and correlation length exponents are derived. The 
critical exponents satisfy known scaling relations and agree well with analytical predictions. 
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\ Driven elastic manifolds in disordered media model a multitude of physical systems ranging from cracks advancing 
in a solid in response to an external strain field [l| to charge density waves |2|, vortices in type-II superconductors 
in an electric field 0], and to magnetic domain walls in an external magnetic field These widely studied out-of- 
equilibrium systems undergo dynamic phase transitions as a function of the external driving force (strain, electric 



or magnetic field, etc.). At driving forces below the depinning threshold (/ < / c ), the disorder pins the manifold so 
that its center-of-mass velocity vanishes at large times. Above the depinning threshold f c , avalanches advance the 
manifold with finite center-of-mass velocity. At the depinning threshold (for / — * /+), the center-of-mass velocity 
in the long-time limit goes to zero algebraically, whereas the typical length, width and duration of the avalanches 
diverge. This critical divergence is characterized by universal scaling exponents. 

Most simply, advancing crack fronts are modeled by the dynamics of an elastic string driven through quenched 
\ disorder (the inhomogeneous material) . On that level, elasticity theory combined with a reasonable energy dissipation 
mechanism can be reduced to an effective equation of motion for the crack front [l], 0, 0] ■ The resulting elastic energy 
kernel is long ranged, i.e. falls of as the squared inverse distance. The same long-range elastic energy kernel is found 
in the wetting problem of an advancing contact line Q . 

In the last fifteen years, well-controlled experiments have been designed to measure universal scaling exponents 
HH [I2I [l3| . But the experimental precision available does not yet allow to give definite answers on questions 
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about universality. In addition, experimental results on the roughness exponent of contact lines are inconsistent with 
theoretical predictions. However, a l ong -standing controversy on numerical values of the roughness exponent has 
(T) ! recently been resolved 0, [U, H [U E& (3 E3| ■ 



In this paper, we perform a numerical simulation of the depinning phase transition of the long-range string driven 
through quenched disorder, with the aim of computing the critical exponents to high precision. Numerical simulations 
of this model are non-trivial because of large finite-size effects, caused by the unusually large value of the finite-size- 
scaling exponent. 

The discrete long-range elastic string is represented by a vector of integer heights {hi, h^} with periodic bound- 
ary conditions. The string follows discrete quasi-static dynamics while being driven through a discrete disorder (cf 
[2l|). The long-range elastic force acting on a string element i is given by 



^elastic n j ' "» (1 \ 



Even though the sum contains ~ L terms for each of L sites, all forces on all sites can be computed in ~ L log L 
operations using a fast-Fourier-transform (FFT) algorithm. 
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The quenched disorder generates a random force r\ 



fr d ° m = f + 9rn(h). (2) 

While / is a constant external field, the random position-dependent force 77 is drawn from a binomial distribution: 
•q equals plus one with probability p, and minus one with probability 1 — p. The coupling g [22I ] controls the Larkin 
length at which the elastic and disorder forces balance. 

The random force results in a net external driving force of / + g(2p ~ I). The parameter / can be adjusted to keep 
constant the otherwise fluctuating center-of-mass velocity v of the string. But we use fluctuating- velocity /fixed- force 
simulations throughout this paper, unless stated otherwise, p is the control parameter of the dynamic phase transition: 
its distance to the critical value p c (the depinning threshold) controls the large-scale properties of the dynamics, the 
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FIG. 1: At the depinning threshold p c , the ratio T(t) = w/[tv(t)] is constant in time, as shown by scaling arguments (for 
p < p c , Y(t —* oo) — > oo because the velocity vanishes, and for p > p c , Y(t — » oo) — > 0). The inset shows fluctuations around 
zero of the external force / for a simulation with fixed p c and fixed velocity: one percent of all sites move at each time step. 



divergence of the correlation length and the algebraic decay of the mean velocity v [23|. The skewed binomial 
distribution of random forces corresponds to random-field frozen disorder, as opposed to random-bond disorder. 
The microscopic dynamics of the string is given by the following discrete dynamical rule: 



hi(t + I) - hi(t) =Vi(t) 



•£ yelastic _|_ y random ^ q 

otherwise 



f = 0,1,. 



(3) 



This rule does not allow for backward motion, as is justified by the no-return theorem 24 1. 

Very large system sizes are crucial in order to calculate the critical exponents with sufficient precision. Consider for 
example the velocity exponent /3, defined as v ~ (p—p c )P. For it to have an accuracy of 5/3 = 0.01 requires a resolution 
of Sp c = 10~ 3 in the critical force (2||. This resolution can only be obtained if the finite-sample-size fluctuations of the 
critical force <t Pc are small enough. From a finite-size-scaling analysis one finds a Pe (L) = 0.65 x L' 1 ^ 3 , v FS = 1.62 
(data not shown), and demanding that a Pa < 10 -3 gives the minimum system size L m - ln w 10 5 . 

The discretized model, with its coarse and low-cost description of microscopic dynamics, can be simulated for sizes 
up to L = 2 20 ss 10 6 , more than three orders of magnitude larger than for earlier discrete 17, [l8[ and continuum 
simulations (26[. 

To benefit from the high speed of the FFT algorithm, we use parallel updating: The elastic force on all string 
elements is calculated. Then, the dynamic rule (equation ([3])) is applied simultaneously to all sites. Finally, the 
random force is updated for those sites that have stepped forward. In this way a single vector of length L suffices to 
store the random force, avoiding storage problems even for very large system sizes. The parallel updating leads to a 
fluctuating mean velocity v, whereas the external driving force is kept constant via the fixed control parameter p. 

The depinning threshold p c is determined by adapting a technique known from absorbing-state dynamic phase 
transitions (2?]]: we start the simulation for fixed p, f = from a flat line at t — 0, and monitor the ratio T of string 
width w(t) = [(1/L) ^2f{hi — (/i)) 2 ] 1//2 to decaying center-of-mass velocity v(t) = (1/L) Y^f Vi(t) times time t: 



At threshold 



v(t) 



w(t) 



„. . w(t) 

Tit) = — -— ~ constant. 

^ ; tv(t) 



To understand why the velocity at threshold is given by the width divided by time, we remark that, intuitively, an 
avalanche takes time t to advance a distance equal to its width w(t). Moreover, at threshold, the motion of the string 
is solely due to avalanches, hence the mean velocity is given by v(t) = w{t)/t. It follows that the ratio T(i) is constant 
at threshold. Away from threshold it either grows to infinity (for p < p c ) or decays to zero (for p > p c ), providing a 
very sensitive indicator for the depinning threshold. 

The main plot in figure [1] illustrates the sensitivity of T(i) to tiny changes in p — p c , which in turn enables a very 
precise localization of the critical threshold p c — 0.5931(2). p c can alternatively be extracted from a power-law fit of 
the steady-state mean velocity (v ~ (p — p c )P , see figure [5]). 
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FIG. 2: Roughness exponent At short distances, the two-point correlation function scales as C(x) oc x 2 ^ . A linear fit gives 
£ = 0.385(5) (average taken over pinned configurations exactly at p — p c ). The inset compares correlation-function data (for 
L — 1024) with the Gaussian model of (eq. [4]). 



There has been debate about the equivalence between simulations with fluctuating velocity /fixed force on the one 
hand and with fixed velocity/fluctuating force on the other. In the thermodynamic limit, one expects them to be 
equal because the fluctuations of the center-of-mass velocity vanish on large enough length scales. In experimental 
or numerical systems, the center-of-mass velocity fluctuates only on scales smaller than the correlation length. In 
our very large simulations, we are able to check the equivalence between the two approaches: in addition to the 
fluctuating-velocity simulations in the rest of this paper, we perform a fixed-velocity simulation by adjusting the 
parameter / at each time step such that a very small but fixed proportion of sites move, thus automatically tuning 
the system to the dcpinning threshold. The inset of figure Q] shows the minute fluctuations of the external force / 
around a zero mean value in a simulation where the external force is chosen at each time-step such that one percent 
of sites are made to move (28|, at p = p c . This demonstrates that p c is indeed the critical point, and that the two 
types of simulations are equivalent. 



The roughness exponent ( = 0.388 |20j characterizes the self-affine spatial structure of the string at depinning, and the 
scaling properties of its average width w (see also j29] ) . As a check for the discrete model, we compute this exponent 
from the short-distance behaviour of the two-point correlation function C(x) = ((h(0) — h(x)) 2 ) ~ x 2 ^ . Averaging 
over pinned configurations exactly at the depinning threshold p = p c , we confirm earlier results (£ = 0.385(5), see 
figure . 

The correlation function C(x) can be fitted over the whole interval x S [0,0. 5L] using the Gaussian model [30| : 
suppose the random string configuration h(x) is decomposed into Fourier modes 



h{x) 



a n cos(2imx / L) + b n sm{2imx/ L), 



where a n and b n are Gaussian random variables with variance a 2 oc l/n 1+2t '. The resulting two-point correlation 
function is given by C Gauss (x) = n sia 2 (imx/L)/n 2< ' +1 , with the following expansion for small x/L [3(| (note the 
difference between the roughness exponent £ and the Riemann zcta- function £r): 



C Gauss (a;) = constant j - (|) 2? 2- 1+2C 7r 2C r(-2C) sin 



tt(2C + 1) 



x\ 2 2tt 2 
LJ ~2T 



- 1) - (f) 4 ^C- 3, + [(f)«]} 



(4) 



This function is used in the inset of figure [2] to fit the correlation function data obtained with the exact algorithm 
3l| for smaller systems. The fit remains excellent even for large x/L. 
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FIG. 3: Growth exponent £/z. The algebraic growth of the average width w(t) ~ ^ , at the depinning threshold (p = p c ) and 
after initial relaxation effects (t ~ 100), gives the growth exponent C,/ z = 0.495(5). Inset: large finite-size effects as width levels 



off at wl ~ L c , for L = {2 1 



,2 20 }. 



The dynamic exponent z was previously estimated for discrete extremal models of limited size [171 Il8| . We determine 
it via the growth exponent (/z, which controls the growth of the average width with time w(t) ~ t^l* . This algebraic 
growth law results from the two scaling laws for the width and the time in terms of the correlation length £: w ~ $ 
and £ ~ t x l z . We start at t = from a flat line, at the critical point p = p c , and monitor the growth of the average 
width. Figure [3] shows how, after an initial relaxation phase [32j|, the width grows algebraically with the growth 
exponent (/z = 0.495(5), implying z = 0.770(5) for the dynamic exponent. The width w(t) is also affected by large 
finite-size effects. The algebraic growth regime is limited by the width of the entire string wl ~ L v The inset of 
figure [3] demonstrates how w(t) saturates at the system-size-dependent limit wl- The plots for different system sizes 
L can be superposed by plotting w(t)L~^ versus tL~ z (not shown). 

Another way of determining the dynamic exponent z is to analyze the time-dependent structure factor S(q,t) = 
(h(q,t)h(—q,t)) with the Fourier transform of the height given by h(q) = X)n=o exp(ign)/i„ and q = 2irk/L, where 
k G {0, . . . , L — 1} is the mode number. The structure factor allows us to illustrate both the dynamic growth of the 



interface and its self-affine structure in space-time 33 1 



f 

S{q,t) ~ q^^t) with *(„) = rV'f • (5) 

I constant , y — + cxd 

Figure 2] illustrates the convergence of the different Fourier modes S(q,t) to their steady-state value ~ q^ 1 ^ 2 ^ after 
a relaxation time <7 -z . The small-time behavior S(q,t) ~ f( 1 + 2 f)/ 2 provides the dynamic exponent, and we confirm 
the numerical value z = 0.77(1), albeit with lower precision than in our analysis using w(t). The inset of figure 0] 
plots S(q,t)q 1+21 ^ versus the dimensionless quantity q z t. It shows the universal scaling form <f>, and demonstrates the 
self-affine link between time t and (inverse) space q, via the dynamic exponent z. 

The velocity exponent (3 controls the power-law decrease of the mean velocity with the distance to the critical 
threshold: v ~ (p — Pc) 13 ■ It is obtained from the center-of-mass velocity v in the steady-state regime (see figured]). 
From a fit at fixed p c , we obtain f3 = 0.625(5). The data of the steady-state velocity can also be fitted to obtain the 
threshold p c , as shown in the inset of figure[5l Only for p = p c is there a straight line on the log — log plot. For p ^ p c , 
a clear non-zero curvature emerges. 

Note that it is quite difficult to obtain a precise estimate for /3, again because of large finite-size effects. Trying 
to calculate (3 on medium-size systems one faces the problem of a significantly reduced critical scaling regime. It 
manifests itself already in the short-range elastic model 26], and is even more pronounced on long-range elastic 
interfaces of intermediate size (data not shown). Our very large systems allow to obtain satisfying accuracy. 
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FIG. 4: Time- dependent structure factor S(q,t). Different Fourier components [from (+) q = 3 x 10 -6 to (•) q = 10 -2 ] take 
a relaxation time of q~ z to reach their steady-state value of q ,_1_2l >. Inset: the small-time behaviour of the universal scaling 
function gives a further estimate of the dynamic exponent z = 0.77(1). 



Our value for the velocity exponent yields an independent check of the growth exponent (,/z, via two sca ling 



relations that reduce the set of four critical exponents to two independent ones: the statistical tilt symmetry 34 1 
links the roughness exponent to the correlation length exponent v = 1/(1 — C), and a further scaling relation links all 
four exponents (3 = v(z — £) (HI. Using these two equalities and values of C = 0.385(5) and /3 = 0.625(5) we converge 
on a value for the growth exponent of Q/z = 0.500(5). These numerical values, our final result, reflect our numerical 
measurements and the constraints of the scaling relations, which are finally used to deduce the dynamic exponent 
z = 0.770(5), and the correlation length exponent v = 1.625(10). 

In conclusion, we have numerically studied the non-equilibrium depinning phase transition of the long-range elastic 
string driven through quenched disorder. The values obtained for the critical exponents satisfy scaling relations and 
agree with analytical predictions from functional RG 13, HI, 3|- For completeness we recall the one- and two- loop 
(e,e 2 ) results: Ce = 0.33, = 0.47; /3 £ = 0.78, (3 e 2 = 0.59; z e = 0.78, z e 2 = 0.66; v e = 1.33, ^ = 1.58. On the 
experimental side, it will be interesting to understand why exactly the roughness and growth exponents calculated 
in this paper agree with those measured in a recent crack front experiment (35| . whether they can also be seen in 
wetting experiments, and more generally, whether universality can be confirmed as the physical mechanism unifying 
these diverse phenomena. 

We thank D. Bonamy, E. Bouchaud, H. Chatc, A. Kolton, P. Le Doussal, S. Guibert, L. Ponson, A. Rosso, and 
K. J. Wiese for stimulating discussions. 




FIG. 5: Velocity exponent f3. The steady-state mean velocity as a function of the control parameter v ~ (p — 
velocity exponent: ft = 0.625(5). Inset: slightly wrong values of p c lead to non-zero curvature ((x) p c — > p c — 

Pc + Spc). 
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